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ABSTRACT 


Appropriate  thermodynamic  models  and  thermo-chemical  data  for  multicompo¬ 
nents  and  immiscible  phases  have  been  Incorporated  Into  a  code  for  the  compu¬ 
tation  of  chemical  equilibrium  of  the  reactants  used  In  a  stored  chemical 
energy  propulsion  system.  The  reactants  considered  are  molten  lithium  fuel  and 
gaseous  sulfur  hexafluoride  oxidant.  Extensive  equilibrium  state  relation¬ 
ships  as  a  function  of  mixture  fraction  have  been  presented  over  a  temperature 
range  of  possible  operating  conditions.  They  were  subsequently  employed  to 
predict  the  complex  flow  structure  of  a  choked,  turbulent,  reacting  SF^  gas 
Jet  submerged  In  a  molten  lithium  bath.  Detailed  profiles  of  velocity,  tem¬ 
perature,  void  fraction  and  mass  fraction  of  all  phases  In  the  Jet  have  been 
presented..  The  results  show  that  the  reaction  Is  completed  within  a  short 
dlstancefrom  the  Injector  but  the  plume  Jet  penetrates  farther  beyond  due  to 
significant  evaporation  of  lithium  fuel  In  the  hot  temperature  plume,  and 
subsequent  condensation  of  the  evaporated  fuel  and  reaction  products  In  the 
/later  part  of  the  plume.  The  effect  of  lithium  subcooling  on  flow  structure 
^  has  also  been  examined.  It  Is  found  that  lowering  the  degree  of  subcooling  by 
raising  the  fuel  bath  temperature  Increases  fuel  evaporation,  decreases  the 
entrainment  rate,  delays  the  completion  of  reaction  and  lengthens  greatly  the 
Jat  penetration  distance.  Finally,  for  practical  use,  a  simple  correlation  of 
the  Jet  penetration  length  as  a  function  of  subcooling  has  also  been  presented. 
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NOMENCLATURE 


Interpreted  as  Interaction  energy 
quasi-lattice  parameters 
specific  heat 
Injector  diameter 
mixture  fraction 

the  square  of  mixture  fraction  fluctuation 
heat  of  reaction 
turbulent  kinetic  energy 
plume  length 

moles  of  chemical  species  1 

gas  constant 

Temperature 

axial  velocity 

axial  distance 

mass  fraction 


Greek  Symbols 

Y|j  activity  coefficient  of  kth  species 

p  density 

c  turbulent  dissipation  rate 


Superscripts 

~  Favre  averaged  quantities 

time  averaged  quantities 


x  - 


Subscripts 

k 

0 

s 

L 

G 

d 

u 


k  species 

jet  nozzle  condition 

bath  (or  surrounding)  condition 

solid 

liquid 

gas 

down  stream 
up  strea.n 
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1.  INTRODUCTION 


Afttr  an  extended  period  of  study.1'5  the  closed  Ranklne-cycle,  lithium 
fueled  sceps  (Stored  Chemical  Energy  Propulsion  System)  his  emerged  es  the 
best  choice  for  undersea  propulsion  system.®  ^  The  SCEPS  derives  Its  eheml- 
cal  heat  source  from  the  highly  exothermic  reaction  between  an  alkali  liquid 
metal  fuel  and  a  halogen  gas  oxidant.  The  reactants  are  so  selected  to  yield 
several  advantages  and  unique  features  for  underwater  propulsion,  first,  the 
energy  density  value  Is  exceptionally  high  on  a  total  reactants  (fuel  plus 
oxidant)  gravimetric  basis.  Secondly,  the  products  of  combustion  are  In  con 
densed  (liquid  and  solid)  phases  and  their  combined  density  Is  greater  than 
the  density  of  molten  fuel.  These  characteristics  permit  the  use  of  e  closed 
reactor  for  a  near-constant  volume,  bath-type  combustion  process,  further¬ 
more,  the  reactor  initially  holding  the  molten  metal  fuel  will  become  the 
holding  tank  to  house  combustion  produett  and  no  discharge  of  products  from 
the  reactor  to  the  environment  Is  necessary,  Coupling  this  thermal  energy 
source  reactor  to  a  closed  Ranklne  steam  cycle  to  propel  a  turbine,  the  per¬ 
formance  of  such  a  combined  power  system  Is  therefore  Insensitive  to  operating 
depth,  making  It  well  suited  for  undersea  app  Meat  tons . 

Among  reactant:  for  undersea  propulsion,  the  use  «f  lithium  (It)  fuel 
with  sulfur  hexafluoride  (SF^)  oxidant  appears  to  be  the  best  candidate  for 
the  SCEPS  concept.  They  ere  relatively  economical.  The  lfft  oxidant  is 
Inert  end  non-toxic.  The  heat  of  reaction  based  on  a  unit  mats  of  total  re 
actants  exceeds  even  that  of  H?/02  reactants.1  They  ere  being  adopted 
for  the  on-going  SCEPS  project.  However,  to  the  authors'  best  knowledge,  lh* 
flow  structure  analysis  of  turbulent  SF^  gas  Jets  submerged  In  end  reeling 
with  a  molten  lithium  has  not  been  reported  In  the  literature. 

•  1  - 


Prior  work  on  the  related  subject  Is  very  limited.  A  Na/Cl  system  was 

8  9 

analyzed  by  a  simplified  Integral  model  which  was  Improved  by  the  use  of 
the  local  homogeneous  flow  approximation,  equilibrium  combustion  model  and 
(c-c-g  model. ^  The  predictions  seem  to  agree  with  the  data  of  a  highly 

O 

diluted,  non-evaporating  sodium  fuel.  Later  it  was  incorporated  into  the 

model  a  general  scheme  for  complex  equilibrium  computation  for  nonideal  multi- 

11  l  ? 

phase  electrolytic  system,  which  was  tested  successfully  with  the  data 

of  a  reacting,  evaporating  and  condensing  system— HCl/aqueous  ammonia  solution 

containing  up  to  25%  of  NH_.  .  to  simulate  the  metal  fuel. 

3(aq) 

9  1 1 

This  study  draws  from  the  prior  works  ’  and  extends  to  Li/SF,  sys- 

0 

tern.  Specifically,  the  thermodynamic  properties  of  the  multiphase,  nonideal 

reacting  mixture  of  Li.,./SF.,  .  reactants  in  chemical  equilibrium  are 

(L)  6(g) 

constructed  and  presented  as  a  function  of  mixture  fraction  (state  relation¬ 
ships).  They  are  then  used  to  predict  the  turbulent  flame  structure  of  Li/SF, 

6 

reacting  jets.  The  former  provides  the  fundamental  information  useful  for 
combustion  studies  in  general  while  the  latter  predicts  the  detailed  flow 
structure  of  the  complex  flame-like  jets  which  may  contain  several  liquid  and 
solid  phases  due  to  partial  miscibility  of  molten  lithium  fuel  and  its  reac¬ 
tion  salt  product,  lithium  fluoride. 


2.  EQUILIBRIUM  STATE  RELATIONSHIPS 


2.1  Formulation 

The  information  on  equilibrium  state  relationships  are  prerequisite  to  the 

analysis  of  the  flame  structure  of  a  SF,  gas  jet  reacting  in  a  molten  Li  fuel. 

0 

The  required  relationships  include  equilibrium  compositions  of  all  species  in 

all  phases,  the  equilibrium  temperature  and  density  of  the  mixture  resulting 

from  the  reaction  of  molten  Li  with  gaseous  SF,  at  various  mixture  fractions. 

0 


2 


The  general  formulation  for  the  state  relationships  has  been  described 
earlier.9,11  Briefly,  It  utilizes  a  variation  of  the  direct  Gibbs  free 
energy  minimization  approach.  The  minimization  Is  subject  to  element  abun¬ 
dance  and  charge  balance  constraints.  A  general  purpose  code  was  developed 

for  complex  equilibrium  computation  for  nonideal  multiphase  systems  (CEC-NMS), 

11  13 

which  has  been  extensively  tested  '  for  a  simulant  reactant  pair  of  HC1 
as  oxidant  and  aqueous  ammonia  solution  as  fuel  to  confirm  the  validity  of  the 
overall  approach. 

One  of  the  objectives  In  the  present  study  Is  to  include  In  the  code  the 

specific  thermodynamic  models  and  the  necessary  thermo-chemical  property  data 

of  the  SF  -LI  reacting  system  and  then  perform  the  complex  equilibrium  com- 
0 

putatlon  of  the  SF&  and  LI  reactants  for  SCEPS. 

2.2  Ll-$Fg  Reaction  System 

The  stoichiometric  reaction  between  lithium  and  sulfur  hexafluoride  can 

3 

be  expressed  as 

8  Li  SF  •*  6  L1F  +  L1.S 

0  C 

which  liberates  a  large  amount  of  thermal  energy  and  produces  lithium  fluor¬ 
ide,  the  dominant  product,  and  lithium  sulfide,  the  minor  product.  For  gener¬ 
ality,  whether  the  reaction  is  stoichiometric  or  not,  the  gas  phase  of  the 
reacting  mixture  In  the  present  study  allows  for  the  possible  existence  of 
fifteen  species.  They  are  Li,  Li2,  F,  F2,  LiF,  Li2F2,  LljFj,  SF,  SF2,  3F^,  SF6,S2F2,Li2S,  S 

and  S2>  which  are  assumed  to  behave  as  an  Ideal  gas  mixture.  Four  solid  spe¬ 
cies  are  considered  in  separate  phases:  Ll^,  ^(s)’  S($)  an<*  2S(sr 
The  liquid  phase  Is  assumed  to  contain  three  possible  species:  Ll^, 

LI F^ L j  and  LlgS^j.  The  liquid  metal,  Ll^j,  and  the  molten  salt,  LIF^, 
can  form  two  Immiscible  phases  with  an  additional  component,  L1?S^j,  dissolved 
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in  each  phase  to  form  ternary  mixtures.  The  activity  coefficients  of  the  Im¬ 
miscible  liquid  phases  (metal-rich  phase  (L-l)  and  salt-rich  phase  ( L-2) )  are 

14 

predicted  using  the  multicomponent  van  Laar  model.  The  activity  coeffi¬ 
cient  of  the  kth  species  may  be  given  as 


In  Y, 


l  l 
1  J 


(A 


Ik 


*  0-5*1]) 


^N^/IRT  (X  N^n 


(D 


t  I  I 

where  and  b^  are  empirical  parameters.  Note  that  A^  =0,  and  A^  * 

I 

Ajj.  The  equation  given  here  differs  from  that  presented  by  Groff  and 
Faeth  (ref.  15,  p.  22,  and  ref.  16,  p.  328).  We  believe  this  equation  to  be 
the  correct  one. 

For  a  system  with  three  components,  the  number  of  empirical  parameters 
may  be  reduced  by  one  by  defining  new  parameters  A^  and  8^  as1*: 


‘u 


b1Ai j ' 


bl/b2: 


b2/b3; 


and  B, 


Vs 


b,b2 


(2) 


The  subscripts  1,  2.  and  3  refer  to  11^.  lIF^j,  and  Li2S(L)’  resPcct,ve1y* 
Equation  (1)  may  then  be  written  as 


Ln 


Yk  =  X  l  (A1k  -  0.5A1J)  (bj/bk)  N^j/RT  [X  Ni  (b1/bJc> 


The  estimated  empirical  parameters,  A^,  from  experimental  data  in  units  of 


(kj/mole)  are  a  function  of  temperature  (K)  as. 
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A12  =  121.53  -  0.05730  T;  A]3  =  53.287  -  4.00E-4  T;  A23  - 
B1  =  2.5846  -  8.132E-4  T;  B3  =  1.36;  B?  *= 


(3) 


(4) 

(5) 


2.3  Thermodynamic  Oata 

Thermodynamic  data  for  all  three  forms  of  H2S  Is  extremely  limited. 

1 5 

Groff  and  Faeth  estimated  the  required  spectroscopic  data  for  determining 
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the  properties  of  L1.S.  .  using  statistical  methods.  They  felt  the  accur- 

2  (g) 

acy  of  the  resulting  thermodynamic  data  for  Ll^  ^  to  be  within  ±  5%.  Their 
equations  are  used  in  this  work  to  calculate  the  enthalpy  and  Gibbs  energy  of 
gaseous  lithium  sulfide.  They  also  estimated  thermodynamic  data  for  ^^(s) 
and  LljS^j  from  Na2S,  1120,  and  Na20  using  Kopp's  rule.17  Their  estimated 
values  were  used  to  calculate  the  enthalpy  and  Gibbs  energy  for  solid  and  liq¬ 
uid  lithium  sulfide.  The  thermodynamic  data  for  all  other  species  were  taken 

1 8 

from  the  data  file  provided  with  the  CEC-72  computer  program. 

2.4  Comparison  with  Experimental  Oata 

The  equilibrium  results  to  be  presented  later  Involve  a  multiphase  and 
multi-species  system.  No  experimental  data  Is  available  to  compare  these 
results  to,  nor  have  any  calculations  of  this  type  appeared  In  the  literature. 
It  Is  decided  to  test  the  validity  of  the  above  thermodynamic  model  and  the 
accuracy  of  the  numerical  calculations  by  computing  the  phase  equilibrium  dia¬ 
grams  of  binary  systems  of  which  experimental  data  are  available  for  compari¬ 
son.  The  tests  Include  LI-ll^S,  11-llF,  and  LIF-ll^S  binary  systems.  A 
representative  comparison  between  the  prediction  (Indicated  by  dots)  and  ex¬ 
perimental  data1*’19  (line)  of  the  LI -LI 2S  binary  Is  shown  In  Fig.  1. 

Good  agreement  was  found  In  all  binaries. 

2.5  Equilibrium  Results 

The  above  thermodynamic  model  and  data  have  been  properly  Incorporated 

into  CEC-NMS  code  to  calculate  equilibrium  state  relationships.  As  stated 

above,  the  resulting  code  has  been  tested  against  the  binary  data  and  Is  now 

ready  for  multi-species  computations.  The  equilibrium  results  are  presented 

In  Figs.  3  to  6  for  the  reaction  of  a  gaseous  SF.  at  T  -  298  K  with  a 

o  o 
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molten  LI  at  T  *  1130  K,  1358  K,  1415  K,  1586  K,  respectively  and  one 
0 

atmosphere  pressure. 

Figure  3  to  6  show  the  equilibrium  temperature,  density,  void  fraction, 

mass  fractions  of  gas,  liquid  and  solid  phases  as  a  function  of  mixture 

fraction  f.  The  latter  Is  so  defined  that  f  •>  1  designates  the  pure  SFC  gas 

and  f  -  0  Indicates  the  pure  molten  LI  metal.  In  Figs.  7  to  34  the 

corresponding  mass  fractions  of  all  Individual  species  In  each  of  three  phases 

are  provided.  By  referring  to  Figs.  3  to  6,  as  the  value  of  f  In  the  reacting 

system  decrea«es  from  one,  the  equilibrium  temperature  Increases  while  the 

mixture  density  decreases.  However,  the  temperature  Is  still  low  enough  that 

a  condensed  solid  phase  of  LiF^  emerges  in  a  narrow  range  and  another 

condensrl  liquid  phase  of  LIF^  also  appears  In  a  narrow  range,  which  can 

be  verified  by  checking  Figs.  11  to  14.  The  co-existed  gas  phase  Is  composed 

of  the  excess  SF  gas  and  vapors  of  combustion  products.  As  f  decreases  to 
o 

the  stoichiometric  value  of  f$t  -  0.72,  where  temperature  reaches  the 
maximum  and  density  the  minimum,  all  species  are  turned  Into  gaseous/vapor 
phase.  Further  decreasing  in  f  will  result  In  a  lower  temperature,  higher 
density  and  the  appearance  of  a  liquid  phase.  Since  molten  Li  and  L1F  liquids 
are  partially  Immiscible,  which  Is  accounted  for  In  the  present  formulation, 
one  might  expect  co-existence  of  two  liquids:  metal  rich  phase  (L-l)  and 
salt-rich  phase  (L-2).  However  our  equilibrium  results  reveal  only  one  liquid 
In  existence  at  a  time.  This  can  be  explained  as  follows-. 

The  necessary  condition  for  two  liquids  to  co-exist  is  that  the  mixture 
Is  fuel-rich,  l.c.,  when  0  <  f  <  f  .  However  In  most  of  region  In  f  =  0  to 
f$t,  the  mixture  equilibrium  temperature  exceeds  the  consolute  temperature 
{1600  K),  the  maximum  temperature  above  which  two  Immiscible  liquids  can  not 
co-exist.  Naturally,  only  single  liquid  exists.  For  the  rest  of  the  region. 
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LI F  Is  below  the  solubility  limit  so  that  again  only  one  liquid  exists.  A 

16  19 

clear  picture  can  be  seen  by  referring  to  the  Li -LI F  phase  diagram  ' 
shawn  In  Fig.  2.  The  equilibrium  states  In  the  present  study  all  lie  In  the 
one  liquid  region  above  the  two-liquid  dome.  For  this  reason,  the  computer 
printouts  do  not  discern  (L-l)  liquid  from  ( L-2)  liquid  and  similar  species 
are  lumped  together.  For  example,  In  Figs.  23  to  26,  the  liquid  lithium  fuel 
species  Is  then  Indicated  by  (L1^  ^  +  Li(L-2)>  or  ^L1(L-1)^'  and  s0 
forth. 

From  the  results  of  the  current  equilibrium  state  relationships,  one  can 
expect  that  no  two  Immiscible  liquids  will  co-exist  In  the  reacting  plume  to 
be  studied  In  the  later  section.  However,  It  does  not  preclude  the  two  Im¬ 
miscible  liquids  from  forming  outside  the  plume  region  when  the  mixture 
temperature  drops  because  of  heat  losses  to  surrounding  heat  exchangers. 

It  Is  Interesting  to  point  out  that  for  the  Identical  Initial  oxidant 

temperature  of  298  K  and  the  fuel  temperature  of  1130  K,  an  Immiscible  liquid 

g 

region  does  appear  In  sodium-chlorine  reacting  system,  but  not  In  the  pre¬ 
sent  lithium-sulfur  hexafluoride  system.  Had  the  Initial  bath  temperature  of 
lithium  been  lower.  It  Is  likely  that  the  solubility  limit  of  the  salt  In  the 
metal  would  have  been  reached  at  a  lower  temperature.  Then  the  salt-rich  liq¬ 
uid  would  begin  to  form  at  a  temperature  below  the  consolute  temperature  and 
an  Immiscible  liquid  region  would  appear. 

Similar  equilibrium  results  have  been  obtained  for  lithium  fuel  maintained 
at  four  elevated  temperatures  listed  In  Table  1.  They  will  be  used  later  In 

the  study  of  subcooling  effects  on  the  flow  structure  of  SF  jets  reacting 

0 

In  molten  Li  pools. 
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3.  STRUCTURE  OF  SF.  JETS  IN  MOLTEN  Li 

o 


3.1  Mathematical  Model 

The  combustor  under  consideration  Is  a  large  vessel  filled  with  a  molten 
LI  fuel  at  one  atmospheric  pressure.  A  gaseous  oxidant  SF  drawn  from  a 

O 

liquefied  storage  tank.  Is  discharged  Into  a  large  vessel  through  a  small  In¬ 
jector  mounted  at  the  bottom  of  the  vessel.  The  choked,  turbulent  SF  gas 

0 

jet  will  react  with  the  entrained  molten  LI  to  form  a  complex,  multiphase, 
ax  1  symmetric  jet.  A  steady  jet  Is  the  subject  under  Investigation.  In  prac¬ 
tice,  the  fuel  vessel  Is  finite  In  size  and  the  reaction  Is  transient  In  na¬ 
ture.  The  current  steady  state  calculation  provides  results  at  a  given  time. 
With  the  state  relationships  developed  above,  the  structure  of  the  steady 

reacting  jet  Is  to  be  predicted.  The  flow  model  adopted  here  follows  closely 

9  10 

the  method  developed  previously.  *  Basically,  It  employs  the  local  homo¬ 
geneous  two-phase  flow  model,  the  local  chemically  equilibrium  combustion 
model  with  the  foregoing  equilibrium  state  relations,  and  Shvab-Zeldovlch  type 

of  formulation.  The  conserved  equations  to  be  solved  are  summarized  In  the 
9  11 

previous  papers  ’  and  will  not  be  re-iterated  here.  They  include  the 
Favre-averaged  parabolic  form  of  transport  equations  for  continuity,  momentum, 
mixture  fraction  (f),  turbulence  kinetic  energy  (k),  the  rate  of  turbulence 
dissipation  (e)  and  the  square  of  the  mixture  fraction  fluctuation  (g).  The 
mean  scalar  properties  other  than  those  from  the  solu-  tlon  of  transport  equa¬ 
tions  are  determined  locally  In  the  plume  by  Integrating  the  product  of  the 
corresponding  equilibrium  state  relationships  and  a  clipped  Gaussian  probabil¬ 
ity  density  profile  over  the  entire  range  of  the  mixture  fraction. 
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3.2  Comparison  with  Experimental  Data 

Again  as  In  the  state  relationships,  no  experimental  data  nor  numerical 

calculations  are  available  for  comparison  with  the  flow  structures  of  Li-SFfi 

system  to  be  presented  next.  However,  the  above  mathematical  model  and  the 

resulting  modified  GENMIX  code  developed  In  this  laboratory  were  extensively 

1 2 

tested  using  experimental  data  of  metal  simulant  reactants  (HC1-NH, 

11  13 

aqueous  solution)  In  our  recent  studies  *  Good  agreement  was  found  be¬ 
tween  predictions  and  data  for  various  fuel  temperature,  pressure  and  concen¬ 
tration  conditions.  The  primary  purpose  of  those  studies  Is  In  fact  to  lay  a 
ground  work  for  the  present  liquid  metal  combustion  study,  the  results  of 
which  are  given  below. 


3.3  Li -SF&  Flow  Structure 

The  Initial  condition  of  the  sulfur  hexafluoride  gas  jet  at  the  Injector 

11  13 

exit  Is  assumed  to  be  uniform  and  taken  as  ’ 


Pq  •  5.97  Kg/m3,  Tq  -  298  X,  uQ  *  203.4  m/s 

f0  -  1.  kQ  =  (0.03  uQ)2.  cQ  -  (0.03  uo)3/0,  gQ  -  0 


(6) 


The  Injector  diameter  0  =  7.9  x  10  m  Is  used  In  the  computation.  At  the 

plume  center  axis,  the  radial  gradients  of  u,  f,  k,  e,  and  g  are  set  to  zero 

due  to  symmetry.  At  the  outer  boundary  of  the  jet,  g  *  0,  u^  -  0.0005  u  and 

0 

rw  ^  2 

f  *  0.0005  f  are  assumed-,  while  k  and  c  are  described  by  k  “  (0.03  uQ)  ,  and 
~  3 

c  *  (0.03  uQ)  /D  at  the  first  axial  step.  At  later  steps  on  the  outer 
boundary, 


k.  -  c..  Ax/u 


u  *  cd 


cu  -  Cc2  ‘u  Ax/(ku  V 
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which  are  obtained  by  solving  the  governing  equations  with  the  radial  convec¬ 
tive  and  diffusive  terms  neglected.  The  subscripts  d  and  u  denote  the  down¬ 
stream  and  upstream  points,  respectively.  The  numerical  computational  details 
are  given  In  the  previous  study.11 

Using  the  above  Initial  and  boundary  conditions,  and  Incorporating  the 

equilibrium  state  relationships  In  the  modified  GENMIX  code,  the  structure  of 

a  reacting  SF  gas  jet  submerged  In  a  molten  lithium  bath  has  been  computed. 

0 

F1gi,35  -38s hows  the  predicted  mean  quantities  along  the  center  axis  of  the  jet 

when  the  bath  temperature  Is  maintained  separately  at  T  *  1130  K  and  etc.  , 

00 

In  which  1  atm.  pressure  Is  maintained.  Centerline  velocity  (u/uo),  mix¬ 
ture  fraction  (f),  temperature  (?),  gas  mass  fraction  (Yg),  liquid  mass  frac¬ 
tion  (Yl),  solid  mass  fraction  (Y$)  and  void  fraction  (a)  are  plotted 
against  the  dimensionless  distance  (x/0)  from  the  Injector  exit.  It  Is  seen 
that  the  mixture  temperature  rises  sharply  in  a  narrow  range  of  (x/D).  Com¬ 
bustion  products  appear  in  solid  and  liquid  phases  with  some  overlapping  near 
the  Injector  due  to  relatively  low  mixture  temperature.  Then  they  turn  Into 
vapor  phase  around  the  high  temperature  zone  and  finally  Into  liquid  phase  far 
away  from  the  Injector.  The  complete  condensation  of  reaction  products  Is  one 
of  the  unique  and  desirable  features  of  the  SCEPS,  which  results  In  a  rela¬ 
tively  short  plume  length  of  I/O  =  5.6.  The  plume  length  here  Is  arbitrarily 
taken  as  the  length  where  the  void  fraction  reduces  to  0.01.  It  Is  the  total 
length  which  Includes  the  reaction,  evaporation  and  condensation  regions.  The 
flow  structure  computation  Indicates  that  the  reaction  region  ends  at  x/0  = 

3.9  where  the  Injected  SF,  oxidant  Is  consumed  almost  completely  (the  mean 

D 

mass  fraction  of  SF,  and  all  Its  derivatives  like  dissociated  elements  etc. 

O 

reduces,  say,  to  0.001  along  the  jet  axis).  As  the  mixture  temperature 
Increases  in  the  reaction  region  the  evaporation  of  the  molten  Li  entrained 


Into  the  jet  stream  from  the  surrounding  bath  becomes  Increasingly  Important. 
The  excess  evaporated  fuel  and  combustion  products  are  eventually  condensed 
downstream.  This  evaporation-condensation  region  partially  overlaps  with  the 
reaction  region;  Its  length  Is  rather  substantial  and  cannot  be  Ignored. 
Further  checking  of  species  presented  In  the  evaporation-condensation  region 
Indicates  that  the  lithium  fuel  vapor  outweighs  the  reaction  product  vapor. 
This  can  be  explained  if  reference  Is  made  to  the  equilibrium  composition  In 
Figs.  11  to  14.  Analogous  to  the  flow  structure  diagram,  the  corresponding 
reaction  region  In  Figs.  11  to  14  is  approximately  extended  from  f  *  1  to  0.67 
at  which  the  oxidant  SF^  gas  and  its  derivative  components  are  practically 
all  consumed  In  the  reaction.  Similarly  the  evaporation-condensation  region 
lies  between  f  =  0.79  to  0.13.  It  Is  clearly  seen  that  In  some 
evaporation-condensation  region  the  mass  fraction  of  fuel  (11  +  L1?)  vapor, 
indeed  far  exceeds  that  of  the  vapor  product  (11F  +  Li 2F2  +  Ll^). 

3.4  Subcooling  Effect 

To  Investigate  effects  of  subcooling  on  turbulent  structure  of  reacting 
jets,  the  bath  fuel  temperature  Is  varied  from  1130  K  to  1586  K  as  shown  in 
Table  1.  It  covers  practically  all  the  possible  operation  range  of  SCEPS. 
Since  the  working  temperature  of  the  reactor  must  be  so  chosen  that  the  reac¬ 
tion  products  are  in  liquid  state,  the  lowest  possible  operation  temperature 
is  limited  by  the  melting  point  of  the  mixture  of  lithium  fluoride  and  lithium 
sulfide,  which  Is  1083  K.  Similarly,  *.nt  maximum  temperature  Is  bounded  by 
the  saturation  (boiling)  temperature  of  lithium  fuel,  which  Is  1609  K. 

Flow  structures  of  all  subcooling  cases  listed  in  Table  1  have  been 
computed.  Figures  35  to  38  show  the  effect  of  subcooling  on  the  plume 
penetration.  When  the  bath  temperature  Is  Increased  from  1130  to  1586  K,  the 
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oling  temperature  (l.e.  saturation  temperature  minus  bath  temperature)  is 

rt-uuced  from  479  K  to  23  K  (minimum  subcooling),  and  the  plume  penetration 

length  is  also  increased  in  an  order  of  magnitude.  In  fact  the  reduction  in 

fuel  subcooling  delays  considerably  the  appearance  of  solid  and  liquid 

phases.  Also  delayed  is  the  temperature  peak  which  implies  that  the  primary 

reaction  occurs  at  further  downstream  from  the  injector.  The  latter  is 

attributed  to  the  reduction  in  the  fuel  entrainment  rate  in  spite  of 

increasing  fuel  evaporation.  Since  Li  vapor  is  lighter  than  the  Li  liquid  in 

the  bath  and  SF,  gas  in  the  main  jet  stream,  the  jet  mixture  density 
0 

decreases  near  the  boundary  as  more  Li  evaporates  from  the  jet  interphase 
boundary.  Consequently  the  turbulent  viscosity  and  the  entrainment  rate  of 
the  ambient  fuel  are  reduced  accord-  ingly.  Similar  phenomenon  has  been 
observed  in  the  simulant  reactant  pair  of  hydrogen  chloride  gas  reacting  with 
ammonia  solution.11 


3.5  Plume  Length  Correlation 

The  total  penetration  length  is  one  of  the  key  parameters  needed  for  the 
design  of  the  liquid  metal  combustor.  It  is  desirable  if  the  plume  lengths 
from  the  above  complex  computation  can  be  correlated  by  a  simple  equation  for 
practical  use. 

With  the  consideration  of  a  simple  entrainment  expression  and  the  energy 
balance  between  the  reaction  heat  release  and  heating  and  vaporization  of  the 
entrained  fuel,  the  plume  length  could  be  grossly  estimated  as12 


AHr  ,Vl/2 

cP  aT  V 


(8) 


where  c  is  a  correlation  constant,  a H  is  the  heat  of  reaction  per  unit 

r 

mass  of  oxidant,  cp  is  the  specific  heat  of  the  bath  liquid  and  AT  is  the 
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degree  of  subcooling  of  the  bath  liquid.  Table  1  lists  most  of  the  relevant 
parameters  including  the  computed  plume  lengths  under  various  subcooling  con¬ 
ditions.  The  L  in  the  above  equation  represents  the  length  over  which  all  the 
reaction  heat  is  dissipated  by  the  sensible  heat  of  the  entrained  liquid  from 
the  bath.  It  is  consistent  with  the  notion  that  the  plume  mostly  consists  of 
fuel  vapor  produced  by  reaction  heat  release  and  that  the  plume  length  is 

determined  by  the  distance  required  for  the  condensation  of  the  fuel  vapor. 

20 

The  equation  is  essentially  identical  to  the  correlation  of  Weimer  et  al. 

for  their  data  on  penetration  lengths  of  condensing  vapor  jets  in  liquids. 

20 

The  c  value  proposed  by  Weimer  et  al.  Is  17.8  which  was  also  used  for 
8  12 

Na-Cl2  system.  ’  For  HC1-NH3  aqueous  solution  system,  about  half  of 

12 

that  value  was  suggested.  In  the  present  Li-SFfi  reacting  system,  we 
attempt  to  correlate  the  subcooling  effect  by  equation  (8).  The  value  of  the 
correlation  constant  is  found  to  be  5.506  by  the  least  square  method.  The 
absolute  deviation  in  the  plume  length  of  the  correlation  from  the  computed 
length  is  between  1.8  to  57%. 


4.  CONCLUSIONS 

The  liquid  metal  combustion  in  stored  chemical  energy  propulsion  system 
has  been  analyzed.  The  necessary  thermodynamic  data  have  been  presented  and 
incorporated  into  a  general  purpose  code  to  account  for  nonideal  solutions 
forming  immiscible  phases.  Specifically,  equilibrium  state  relationships  of 
the  gaseous  SF^molten  Li'  reacting  system  have  been  presented  as  a  function 
of  mixture  fraction.  Comparisons  have  been  made  with  limited  experimental 
data  available  and  good  agreement  is  found.  Using  the  state  relationships  and 
appropriate  flow  models,  the  structure  of  sonic,  'turbulent,  reacting  SF  gas 
jets  submerged  in  molten  Li  baths  has  been  obtained.  Velocity,  temperature, 
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void  fraction,  and  mass  fractions  of  gas,  liquid  and  solid  phases,  etc.  in  the 
reacting  jets  have  been  presented.  The  plume  length  is  found  to  be  finite, 
due  to  condensation  of  reaction  products,  but  penetrates  far  beyond  the  reac¬ 
tion  length  because  the  lithium  fuel  evaporated  in  the  high  temperature  region 
must  be  re-condensed  in  lower  temperature  region.  The  subcooling  effect  has 
also  been  examined  and  found  to  greatly  affect  the  flow'  structure  by  shifting 
the  structure  downstream  when  the  degree  of  subcooling  is  reduced.  Finally  a 
simple  equation  is  presented  which  correlates  the  plume  length  as  a  function 
of  subcooling  satisfactorily.-  The  results  of  the  present  study  make  possible 
the  prediction  of  flow  structure  in  the  current  design  of  SCEPS  for  undersea 
propulsion  systems. 
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Table  1.  Subcooling  effects  on  plume  length 


T 

00 

K 

Subcooling 

AT.K 

p» 

Kg/m3 

p0 

Kg/m3 

C 

P 

KJ 

Kg  of  Li 

AH 

r 

KJ 

Kg  of  SF6 

L 

d 

1130 

479 

449 

5.97 

4.13 

19,267 

5.66 

1358 

251 

426 

5.97 

4.11 

19,116 

11.64 

1415 

194 

421 

5.97 

4.10 

19,079 

36.44 

1586 

23 

403 

5.97 

4.10 

18,966 

132.42 
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Fig.l  Solubility  Comparison  Between  Prediction(dots)  and  Data  of  Li  ->  LijS  System 
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Fig. 2  Equilibrium  States  (of  Molten  Li  fuel  at  700=1130  K  reacting  with 
gaseous  SFe  at  70=298  K)  Phase  Diagram 
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Fig. 3  Equations  of  States  For  Li  —  SF6  System  with  Li  at  1130  K  and  5F«  Gas  at  298  K 
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Fig. 4  Equations  of  States  For  Li  —  SF&  System  w:th  Li  at  1358  K  and  SF&  Gas  at  298  K 
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Fig. 5  Equations  of  States  For  Li-  SF6  System  with  Li  at  1415  K  and  SF6  Gas  at  298  K 
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Fig.6  Equations  of  States  For  Li  -  SF6  System  with  Li  at  1586  K  and  SFe  Gas  at  298  k 
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Fig. 23  Equilibrium  Composite 
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Fig. 25  Equilibrium  Compositions  of  Li  —  SFe  System  with  Li  at  1415  K  and  SFq  Gas 
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tem  with  Li  at  1586  K  and  SF$ 
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Equilibrium  Compositions  of  Li  —  SFe  System  with  Li 


Fig.33  Equilibrium  Compositions 
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Fig. 35  Predicted  Mean  Quantities  Along  The  Jet  Axis 
SE 9  Gas  (298K)  Discharged  In  Molten  Li  (1130K) 
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Fig. 36  Predicted  Mean  Quantities  Along  The  Jet  Axis 
SFg  Gas  (298K)  Discharged  In  Molten  Lx  (1358K) 
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MEAN  QUANTITIES 


Fig. 37  Predicted  Mean  Quantities  Along  The  Jet  Axis 
SFt  Gas  (298K)  Discharged  In  Molten  Li  (1415K) 
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MEAN  QUANTITIES 
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Fig. 38  Predicted  Mean  Quantities  Along  The  Jet  Axis 
SF6  Gas  (298K)  Discharged  In  Molten  Li  (1586K) 
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APPENDIX  A.  Equilibrium  Computer  Program  And  Results  For  Li  -  SF*  Reaction 
APPENDIX  B.  Flow  Structure  Computer  Program  And  Result*  For  Li  -  SF^  Reaction 

Legible  copy  can  not  be  reproduced. 

Interested  readers  should  contact  the  author  (S.H.Chan). 
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